Vortex flow around a circular cylinder above a plane 
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The study of vortex flows in the vicinity of multiple solid obstacles is of considerable theoretical 
interest and practical importance. In particular, the case of flows past a circular cylinder placed 
above a plane wall has attracted a lot of attention recently. In this case, a stationary vortex is 
formed in front of the cylinder, in contradistinction to the usual case without the plane where a 
vortex pair is observed behind the cylinder. In the present work, we apply modern complex analysis 
techniques to obtain the complex potential for the problem of one point-vortex placed in a uniform 
stream past a circular cylinder above a plane. A typical streamline pattern is also shown. 



I. INTRODUCTION 

The formation of vortices in viscous flows past cylin- 
drical structures is a problem of considerable theoretical 
interest and practical relevance for many applications [l[ . 
For example, in the case of a flow past a circular cylinder 
a pair of counter-rotating vortices forms behind the cylin- 
der at small Reynolds numbers. This vortex pair then 
goes unstable at higher Reynolds numbers and evolves 
into a von Karman vortex street. This system was first 
studied analytically by Foppl in 1913 Q. Modelling the 
vortex flow in terms of point vortices in an otherwise 
inviscid and irrotational flow, Foppl was able to find sta- 
tionary configurations for a pair of vortices behind the 
cylinder and analyze their stability. (Part of Foppl's 
original stability analysis was however in error, as it has 
been pointed out by several authors [1, H|.) Because it is 
amenable to analytical treatment, the point- vortex model 
is an important tool to study the basic aspects of vortex 
dynamics in real fluids. In addition, the study of vor- 
tex phenomena directly from the Navier- Stokes equation 
is very costly computationally, which makes the point- 
vortex model even more attractive. 

The problem of vortex flow in the presence of several 
solid boundaries (obstacles) is much more difficult, and 
theoretical studies for such cases are more sparse in the 
literature. One geometry of particular interest is that 
of a uniform flow past a circular cylinder placed above 
a planar solid wall, where vortices can form upstream of 
the cylinder Q . Here we wish to investigate this problem 
in terms of a point-vortex model. More specifically, we 
consider a point vortex of intensity T placed in a uniform 
stream of velocity U past a circular cylinder of radius s, 
whose center lies at a distance d from a plane; see Fig. [T] 
The main objective of the present paper is to derive the 
complex potential for this system. A typical streamline 
pattern of the flow will also be presented. 

Since the fluid domain for the case in study (Fig. 1) 
is not simply connected (nor can it be reduced to a sim- 
ply connected region by symmetry arguments), we need 
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FIG. 1. Point vortex in a uniform stream past a circular 
cylinder above a plane. 



to consider conformal mappings between multiply con- 
nected domains in order to compute the complex po- 
tential for the flow. More specifically, we shall make 
use of a mathematical formalism recently developed by 
Crowdy and Marshall @ , which is based on the so-called 
Schottky-Klein prime function, as discussed next. 



II. METHODOLOGY 

We assume that the fluid is inviscid, irrotational, and 
incompressible, so that the fluid velocity v(x, y) is given 
by the gradient of a potential function: v — V(/>, where 
the velocity potential 4>(x, y) obeys Laplace equation: 
V 2 </> = 0. Our main goal here is to compute the com- 
plex potential w(z) = 4>{x,y) + iip(x,y), where ip is the 
so-called stream function, for the problem illustrated in 
Fig. 1. To this end, let us introduce the conformal map- 
ping z(C) from an annular region, ro < \(\ < 1, in the 
auxiliary complex £-plane onto the fluid domain in the 
complex z-plane, where the unit circle |£| = 1 is mapped 
to the plane boundary and the inner circle |£| = tq is 
mapped to the cylinder. Furthermore, the points £ = i 
and £ = — i are mapped to z — and z = oo, respec- 
tively; see Fig. [2] 

It is not difficult to see that the function that enacts 
the desired mapping is given by 
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FIG. 2. Circular domain in the auxiliary complex £-plane 
(left) and physical domain in the complex z-plane (right). 



The radius ro of the inner circle in the £-plane is related 
to the physical parameters s and d by the following ex- 
pression: 



where w(C, a) is the so-called Schottky-Klein prime func- 
tion, which encodes the geometry of the circular domain. 

For the particular geometry shown in Fig. 2, the 
Schottky-Klein prime function w(£, a) can be expressed 
in a simple form: 



where 



W (C,a) = -^P 



c=n(i 

71=1 



c 



-,r 



(7) 



(8) 



and 



1 



ro = 



(3) 



p( x , y) = (i - X ) (i - y 2n x) (i - y 2 "^ 1 ) 



(9) 



as can be easily verified. 

Before presenting the complex potential for the prob- 
lem above, it is instructive to recall [7|] that the complex 
potential for a point vortex of intensity T located at po- 
sition z = z v (in an unbounded domain) is given by 



w(z) = — log (z - z v ) 
2m 



(4) 



Similarly, for a vortex in the presence of a cylinder of 
radius a, the complex potential can be easily obtained 
by virtue of the Milne-Thomson circle theorem [8( , which 
yields 



w(z) = — log 

27TI 



a 2 /z v 



(5) 



Here the term in the numerator comes from the contri- 
bution to the complex potential from the vortex itself, 
see Eq. whereas the term in the denominator corre- 
sponds to the contribution from the vortex image (inside 
the cylinder) , which is necessary to enforce the boundary 
condition that the cylinder surface be a streamline of the 
flow. 

For the geometry shown in Fig. 1, where in addition to 
the cylinder there is an extra boundary (the plane y = 0), 
the circle theorem is no longer of help, and one has to 
resort to an alternative approach in order to compute the 
contribution from the infinite set of vortex images (both 
inside the cylinder and below the plane). In this case, 
it is more convenient to compute the complex potential 
in the auxiliary £-plane and then transform it back to 
the z-plane. Indeed, it can be shown [6|, [9J that for a 
point vortex of unit intensity, located at position ( = a 
in a circular domain in the £-plane, the complex potential 
w v ((,a) is given by 
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From the definition of the complex potential w v ((,a) 
given in Eqs. ([5])-©! one can verify that it satisfies the 
appropriate boundary conditions in the £-plane, namely, 
that the unit circle and the inner circle are both stream- 
lines of the flow. 

The function P(x, y) given in Eq. © is related to the 
first Jacobi theta function t?i, which appears in the the- 
ory of elliptic functions [ioj . Indeed, it is possible to 
derive an expression for the complex potential w v (£,a) 
entirely within the framework of elliptic functions. The 
advantage of the method based on the Schottky-Klein 
prime function is that it can be rather easily extended to 
two-dimensional vortex flows around an arbitrary num- 
ber of obstacles. (This general problem is however be- 
yond the scope of the present work.) 

Using the complex potential for a vortex of unit inten- 
sity given in Eq. (JSJ) , it is possible to obtain the complex 
potential, wjj{z), due to a uniform stream of velocity U 
past a cylinder above a plane (with no other flow ele- 
ments). In this case, it can be shown [ll[ that the cor- 
responding complex potential, wu{C), in the C-plane is 
given by 



wu(0 = -^nUiVd 2 - s 2 ( 



dw v dw v 
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(10) 



with w v as in Eq. ©. Notice that the derivatives in 
Eq. (p~0|) are evaluated at the point a = — i, which is the 
point that is mapped to infinity by the conformal map 
z(C), i.e., z(—i) = oo. In the next section we shall use 
Eqs. (fTUf and ^ to construct the complex potential for 
the vortex flow illustrated in Fig. 1. 



III. RESULTS AND DISCUSSION 

The complex potential w(z) for the problem of a vortex 
of intensity T placed at position z v in a uniform stream in 
the presence of a cylinder and a plane boundary can be 
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FIG. 3. Streamline pattern for the flow studied. The param- 
eters chosen are T — —10, U = 1, and z v — —1.5 + 0.5i. 



The streamline pattern shown in Fig. [3] displays some 
interesting features. Of particular interest is the forma- 
tion of a recirculation zone around the vortex in front 
of the cylinder. A pattern similar to this was seen in 
the experiments on flows past a circular cylinder above a 
plane reported in Ref. Notice that the recirculation 
zone acts as an "obstacle" to the oncoming flow, deflect- 
ing part of the fluid over the top of the cylinder and 
thus reducing the flow though the gap. Note also that 
far above from the cylinder the streamlines tend to be 
straight lines, meaning that far away from the cylinder 
one recovers the uniform flow imposed by the uniform 
stream, as expected. 



obtained by combining Eq. (|10j) with Eq. ([6]), multiplied 
by T. This yields 

w{z) = wu(C(z)) + Tw v (((z), ((z v )) , (11) 

where C,{z) is the inverse mapping given in Eq. ([2]). As is 
well known Q, the velocity field v — (u, v) can be found 
by simply taking the derivative of the complex potential: 

dw 

u-iv = — . (12) 
dz 

This yields a pair of differential equations for x — u(x,y) 
and y = v{x,y), which can be numerically integrated to 
generate the flow streamlines. It is however more conve- 
nient to take an alternative approach, namely, perform 
a contour plot of the stream function ijj(x,y), which can 
be easily obtained by taking the imaginary part of the 
complex potential (fTTj). Each level set ip(x,y) = c, with 
c being a constant, then yields a streamline of the flow. 
Figure [3] shows the streamline pattern obtained from the 
contour plot just described for the case when U = 1, 
r = —10, and the vortex is located at z v = —1.5 + 0.5i. 



IV. CONCLUSIONS 

We computed the complex potential w(z) for a point 
vortex in a uniform stream past a circular cylinder placed 
above a plane wall. From the imaginary part of w(z), we 
obtained the stream function of the flow whose level sets 
yield the streamlines. In particular, we presented the 
streamline pattern for the case when the vortex is placed 
upstream of the cylinder — a situation that is of interest 
to the experiments on viscous flows past a cylinder above 
a plane recently performed by Lin et al. 5]. The next 
step towards a more complete theoretical understanding 
of the system is to study the vortex dynamics in this 
geometry, by allowing the vortex to move with the flow 
velocity (excluding the vortex own contribution) at the 
vortex position. This analysis is currently in progress. 
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